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Abstract 



Random matrix ensembles defined by a mean-field one-body plus a 
chaos generating random two-body interaction (called embedded ensembles 
of (l+2)-body interactions) predict for wavefunctions, in the chaotic domain, 
an essentially one parameter Gaussian forms for the energy dependence of the 
number of principal components NPC and the localization length Ih (defined 
by information entropy), which are two important measures of chaos in finite 
interacting many particle systems. Numerical embedded ensemble calcula- 
tions and nuclear shell model results, for NPC and Ih, are compared with 
the theory. These analysis clearly point out that for realistic finite interacting 
many particle systems, in the chaotic domain, wavefunction structure is given 
by (l-|-2)-body embedded random matrix ensembles. 
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I. INTRODUCTION 



In the last few years the study of quantum chaos in isolated finite interacting particle 
systems has shifted from spectral statistics to properties of wavefunctions and transition 
strengths (for example, electromagnetic and Gamow- Teller transition strengths in atomic 
nuclei, dipole strengths in atoms etc.). Working in this direction, several research groups 
have recognized recently that the two-body random matrix ensembles and their various ex- 
tended versions form good models for understanding various aspects of chaos in interacting 
particle systems In particular using the so called EG0E(l+2), embedded Gaussian or- 
thogonal ensemble of (l+2)-body interactions defined by a mean-field one-body plus a chaos 
generating random two-body interaction, there are now several studies on the nature of oc- 
cupancies of single particle states, strength functions (or local density of states), information 
entropy, transition strength sums and transition matrix elements of one-body transition op- 
erators, Fock-space localization etc., in the chaotic domain of interacting particle systems 
such as atoms 0, nuclei quantum dots 0, quantum computers and so on. Ref. |jl| 
gives a recent overview of this subject. The focus in the present article is on two important 
measures of localization (in wavefunctions and transition strength distributions): (i) num- 
ber of principal components NPC (or the inverse participation ratio IPR); (ii) localization 
length Ih as defined by the information entropy (S*"-^"). It is well established that the NPC 
in wavefunctions characterizes various layers of chaos in interacting particle systems ^j. In 
addition, for systems such as atomic nuclei, NPC for transition strengths is a measure of 
fluctuations in transition strength sums. Similarly the role of Ih in quantum chaos studies is 
well emphasized by Izrailev and more significantly, using nuclear physics examples |^ it 
is well demonstrated that the wavefunction entropy S*"-^° coincides with the thermodynamic 
entropy for many particle systems with two-body interactions of sufficient strength but only 
in the presence of a mean-field, i.e. in the chaotic domain but with a mean-field - therefore 
the significance of EG0E(l+2). Clearly deriving the predictions of EG0E(l-|-2) for NPC 
and Ih are of considerable importance. This problem was addressed in Refs. P,p!0|. In 



results for NPC in wavefunctions, in the so called Breit-Wigner (BW) domain, are derived. 



On the other hand, in |Ty] results in the so called Gaussian domain (the BW and Gaussian 
domains are defined in Sect. II ahead) are derived for NPC and Ih in transition strength 
distributions with only the final results mentioned for wavefunctions. The purpose of the 



present paper is to give a detailed derivation of the results mentioned in for NPC and 
Ih in wavefunctions and subject them to numerical tests. Now we will give a preview. 

Section II gives some of the basic results for EG0E(l+2). In Section III, formulas for NPC 
and Ih in wavefunctions are derived by exploiting the Gaussian nature and the associated 
properties of strength functions in EG0E(l+2). Numerical tests of the theory are given in 
Section IV. Finally Section V gives concluding remarks. 

II. BASIC RESULTS FOR (l+2)-BODY RANDOM MATRIX ENSEMBLES 

Given m fermions in single particle states, assuming at the outset that the many 
particle spaces are direct product spaces of the single particle states, two-body random 
matrix ensembles (usually called TBRE) are generated by defining the hamiltonian H, 
which is 2-body, to be a random matrix in the 2-particle spaces and then propagating 
it to the ^ ^ dimensional m - particle spaces by using their geometry (direct prod- 
uct structure); often one considers a GOE representation in the 2-particle spaces and 
then the TBRE is called EG0E(2); see for more details. For a EG0E(2), with 
N » m » 2, the normalized state density p{E) = {6{H — E)) takes Gaussian form 
and it is defined by its centroid e = (H) and variance = {{H — e)^). In order to ex- 
plicitly state that the state density is generated by the hamiltonian H, sometimes p{E) 
is denoted as p^{E) and similarly e as en and a as an- Note that the averages ( ) are 
over the m-particle spaces ; in the nuclear physics examples, they are usually over the 
m-particle spaces with fixed angular momentum (J) and isospin (T) which are good quan- 
tum numbers. Just as the state density, given a transition operator O, the normalized 



bivariate strength densities (matrix elements of O weighted by the state densities at the 
initial and final energies) pbiv;o{Ei, Ef) = [(^0^0)Y^ {0^5{H - Ef)06{H - Ei)) take bi- 
variate Gaussian form for EG0E(2) and it is defined by the centroids (e^, e/) and widths 
((Tj, cTf) of its two marginals and the bivariate correlation coefficient which is given by 
(Ct [{H - ef)/af] O [{H - ei)/ai]) / (O^O). Thirdly, the level and strength fiuctuations 
follow GOE. Moreover, with the Gaussian form for the state densities and bivariate Gaus- 
sian form for the strength densities, the strength sums (e \ O'^O \ e'^ = J2e' \{E' \ O \ E)f 
take the form of ratio of two Gaussians, (E \ O'^O \ E\ = (0'^0)paio,g{E)/pg{E) where 



Pow.g{E) = {0^06{H - E)) is defined by its centroid e^to = (O^OH) / {O^O) and vari- 
ance = (O'^OH'^^ I iO'^O^ — e^tc?; Q stands for Gaussian. 

Hamiltonians for realistic interacting particle systems contain a mean-field part (one-body 
part /i(l)) and a two-body residual interaction V(2^ mixing the configurations built out of 
the distribution of particles in the mean-field single particle states; /i(l) is defined by the 
single particle energies e^, i = 1,2,...,N and V{2) is defined by its two-particle matrix 
elements. Then it is more realistic to use EGOE( 1-1-2), the embedded Gaussian orthogonal 
ensemble of random matrices of (l-|-2)-body hamiltonians where {H} = h{l) + X{V(2)}; 
sometimes it is more convenient to use ah{l) + A{y(2)}. Here { } denotes ensemble, A 
and a are free parameters and V{2) in the two particle spaces is a GOE with unit matrix 
elements variance; note that in general h{l) need not be fixed nor V{2) a GOE (in this 
general case, the ensemble is called EE(l-|-2); see for more details). At this stage it is 
important to stress that all the EG0E(2) results mentioned before are indeed applicable 
to EG0E(l-|-2) but only in the domain of chaos. Given {m,N) and the average spacing 
A (generated by h{l)) of the single particle states (without loss of generality one can put 
A = 1) it is possible to find the critical A value Ac such that for A > Ac there is onset of 
chaos (GOE level fiuctuations) in many (m >> 1) particle spaces. In fact Ac is of the order 
of the spacing between m-particle mean-field basis states that are directly coupled by the 
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two-body interaction; see the second and third reference in 0]. For A > Ac, for example, 
it is well established that the transition strength sums in EG0E(l+2) follow the EG0E(2) 
forms; see Fig. Ic ahead. Refs. ||1],|1T| give many numerical examples for this, drawn from 
EG0E(H-2) and more importantly for atomic nuclei in several parts of the periodic table 
(detailed discussion of the nuclear physics examples is given in the last reference of [0). It 
should be mentioned that the Gausian forms of state and transition strength densities are 
used in [|l^] to derive simple formulas for NPC and Ih in transition strength distributions. 

For deriving formulas for NPC and Ih in wavef unctions, most useful quantity is the 
strength function (or local density of states) Fk{E). Given the mean-field basis states \k) 
with energies Ei^ = {k \ H \ k), the eigenstates \E) can be expanded as \E) = J2k \k). 
Then the strength function Fk{E) = {6{H - E)f = J2e' Cf ^ 6{E - E') and therefore 
it gives information about the structure of the eigenfunctions. In order to proceed further, 
let us say that the Ej. energies are generated by a hamiltonian (the structure of is 
discussed ahead). With this, it is easy to identify Fk{E) as a conditional density of the 
bivariate density puv{E,Ek) = {6{H — E)6{Hk — Ek)). Taking degeneracies of E and E]^ 
energies into account we have, 

pUE,Ek) = {6{H-E)6{Hk-Ek)) 



E,I3 
k,a 



{l/d)\C^\' dp^{E) dp^^iE, 



Fk{E)= pu.{E,Ek)/p"^{E, 



(1) 



|Cf I' = pUE,Ek)/ Id p^(E) p"^{E, 



In (1) d stands for the dimensionality of the m particle spaces and \C^f is the average of 



over all the degenerate states. Let us now examine the structure of Hk and pbiv{E, Ek). 



Firstly it should be noted that the two-body interaction V{2) can be decomposed into two 
parts V{2) = V^^^+V so that h{l) + V^^^ generates the energies (diagonal matrix elements 
of H in the m-particle mean-field basis states). With m particles in single particle states, 
there is a U{N) group and with respect to this group V^^^^ contains a scalar part (a 
function of m), an effective (m-dependent) one-body (Hartree-Fock like) part V^M'^ and an 
irreducible 2-body part l/M'^. The 1/™'° + 1/™'^ will add to h{l) giving an effective one- 
body part of H; h{l) h{l) + + ^["l'^ = h. The important point now being that, 
with respect to a U{N) norm, the size of yM'^ is usually very small compared to the size 
of h in the m-particle spaces. With this, H = h + V and then the Hk is nothing but h. 
The piece V = V{2) — generates the widths and other shape parameters of Fk{E). It 
should be added that with respect to the U{N) norm h and V are orthogonal and therefore 
ajf = 0"^ + cr|^. Definition of V^^\ a brief discussion of its U{N) decomposition etc. are 
given in Appendix A. For EG0E(l+2), it is well known that the widths of Fk{E) are in 
general constant; see and Appendix A. The average variance of Fk{E)^s is given simply 
by 

4 = ^V = (d)-' E \{^\H\f3)\' 

where a and (3 are m-particle mean-field basis states indices. The results, (i) the norm of 
the part is negligible and (ii) the widths of the strength functions are nearly constant 
(with little fluctuations) are well verified in a number of examples; see and references 
in |I| for many nuclear physics examples. EG0E(l+2) discussions in the literature tacitly 
assume that h is h{l) and V is V{2) and the same is assumed from now on, i.e H = 
h + XV h{l) + XV{2). In addition to (i) and (ii), it is well verified in a number of 
numerical calculations that: (iii) Fk{E)^s exhibit a transition from BW to Gaussian form 
in the chaotic domain defined by A > Aj?^; usually Ac < Xp^'i see flp!^ for some analytical 
understanding of this result. The results (i), (ii) and (iii) clearly imply that the pbiv{E,Ek) 
is a bivariate Gaussian and this result was first mentioned in [jl5[ . A numerical example for 
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the BW to Gaussian transition in strength functions in EG0E(l+2) is shown in Fig. la. 
In this example Ai?^ ~ 0.2 and it is much larger than Ac = 0.06 obtained via the results 
for the Dyson-Mehta A3 level statistic shown in Fig. lb. Thus there is onset of GOE 
fluctuations much before the Fk{E)^s start taking Gaussian form, i.e. ^Ff, > Ac. Unlike 
the case with strength functions (also transition strength densities; see [0), as mentioned 
before, strength sums start following the EG0E(2) form (i.e. pn^,g{E) / pg{E)) from A = Ac. 
This is demonstrated in Fig. Ic where occupancies {E \ni \ E) as a function of E are 
shown (they correspond to the strength sums generated by single state {{i)) destruction 
operators). As mentioned in the introduction, the nature of NPC (which is the inverse of 
IPR) in wavefunctions in the Ac < A < Ai^^. domain where Fk{E) is of BW form (i.e. in 
the BW domain) was studied in p while the present article is concerned with the A > Ai^^, 
domain (i.e. the Gaussian domain) where Fk{E) is of Gaussian form. 

III. EGOE(l+2) RESULTS FOR NPC AND Ih IN WAVEFUNCTIONS 

For EG0E(l+2), in the chaotic domain with A > Xp^, one has from Sect. II the results: 
(i) Ek are generated by = h{l), therefore the variance of p^^^E^) is cr^; (ii) widths of 
the strength functions are constant and they are generated by V{2), the average variance 
a1 = ay] (iii) Ffc(_E')'s are Gaussian in form; (iv) Fk{E) is a conditional density of the 
bivariate Gaussian pbiy,g{E, E^). The correlation coefficient ( of pbiv:g{E, E^) is given by. 

Note that the centroids of the E and Ek energies are both given by en = (H). In (2) the 
second equality is obtained by using the orthogonality between h{l) and V{2) operators. It 
is immediately seen that the is nothing but the variance of E^s (the centroids of Fk{E)) 
normalized by the state density variance. The puv:g{E , Ek) ., which takes into account the 
fluctuations in the centroids of Fk{E) and assumes that the variances are constant, is used 
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to derive formulas for NPC and Ih in the wavefunctions (methods of taking into account 
variance fluctuations will be discussed ahead) ipE — \E) expanded in the mean-field basis 
defined by the states 0^ = \k). Let us first define NPC and Ih, 

\E) = E^f 1^) 

k 



(NPC)^ = 


. k 




1 


(3) 


Ih{E) = 


exp [{S'^^")e 


/(0.48d) ; 






-E \ci 


2 „ 2 

In . 





In (3) 0.48ci is the GOE value for S^"^°, thus Ih = 1 for GOE. Similarly NPC is ci/3 for 
GOE. 



In terms of the locally renarmalized amplitudes Cj^ = Cl^/\\C^f where the bar denotes 
ensemble average with respect to EG0E(l+2), E |cf f = E \^k\* (l^f 1^) ■ Then the 

k A; ^ ^ 

ensemble averaged (NPC) e is obtained as follows, 



E 1^.^ 

k 



4 EGOE{l+2) 



E (nf) 

k ^ ^ 

3 E (^)' 



{3/d) 



I 



pS'iEk) U{E) 



jdEkpf{Ek) [Fk..g{E)f 



(NPC)^ = {d/3) exp-^^^ 



2 b2 ■ 



(4) 



The E in (4) is the standardized i.e. it is zero centered and normalized to unit width, 
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E = {E — eH)/o'H- In the first step in (4) the fact that EGOE exhibits average fluctua- 
tions separation (with little communication between the two) is used. For example, in the 
normal mode decomposition of the EGOE state denity, it is seen that the long wavelength 
parts generate the smoothed Gaussian density (with corrections) and the short wavelength 



parts the GOE fluctuations with damping of the intermediate ones (see [[18H20|| for detailed 



discussions on this important result). This allows one to carry out 



ensemble average 



independent of the other smoothed (average) term. In the second line the Porter-Thomas 



form of local strength fluctuations is used and then |C|'|^=3, a GOE result. In the third step 
the result in (1) is used. Then, the Gaussian forms, valid in the chaotic domain (A > Ai?^.), 
of all the densities for EG0E(H-2) give the final formula (this result was quoted first in 
without details). Before turning to the formula for the localization length Ih, let us briefly 
discuss about the corrections to (4) due to the fluctuations in the variances of Fk{E)] the 
form with Fk{E) shown explicitly, is written in (4) for this purpose and this form also allows 
one to understand the results in as discussed ahead. 

The correction to NFC due to 6al = — (t| 7^ is obtained by using, for small \Sal\, the 
hermite polynomial expansion which gives |T^, Fk.g{E) Fk.g{E) {1 + C2(i^| ~ 1)} where 
C2 = 5a1/2(j1 and £k = {E — Ek)/\[^. This corrected Fk{E) is used in the integral form 
with Fk{E) in (4). As NFC involves sum over all the \k) states, it is a valid assumption to 
treat (5cr|)'s as random and therefore in [Fk{E)]'' only the terms that are quadratic in (Sal) 
will contribute (see P). Replacing [{Sal)/^] by (5^2)/^= [(rf)-^^^ (^^D^Il^^V^ and 
substituting the corrected Fk{E) for Fk,g{E) in (4), we get 

/ dE, p^^m [F,AE)r (1 + ^-^{£1 - 1)V 

>. 1 

XiE) 



(NFC)j 



X{E) 



3 



:i+c^)^ 



C exp 



e E' 



(5, 



cr 



(5) 



3 (1 + 0(1-2(^) ^2 
1-e 



1 + 

1-C^ 



[1 + 2C 



The 5(T^ corection term in (5) is valid only when the fluctuations in the variances of Fk{E)^s 
are small (this is in general always true). For small ( values, the formula for NPC in (5) 
reduces to the expression given recently, for EG0E(2), by Kaplan and Papenbrock ||21]]; they 
use ideas related to 'scar theory'. In the EG0E(l+2) hamiltonian H = h{l) + XV{2), with 
A ^ oo one obtains EG0E(2) and then it is clear from the definition in (2) that in this limit 

2 r/\/\n-l 

C ~ 0. More precisely, with >> m » 1, ~ ( ^ ) and [((5(t^)/(t|'] ~ 



for {H} = {V{2)}; see appendix A. Therefore for finite N, the correlation coefficient and 
the variance corrections are small but non zero and in the large N limit they are zero giving 



the GOE result as pointed out in [|10|. As we add the mean-field part to the EG0E(2), 
( increases and at the same time the variance correction decreases; see Appendix A. Thus 
the formula (5) with the (5cr^) term is important only for small (. Eq. (4) is accurate for 
reasonably large ( (say for ( > 0.3) as in the examples discussed in ||10[. All these results 
are well tested by the numerical examples in Sect. IV. 
Proceeding exactly as in (4), formula for the localization length as a function of the ex- 



citation energy is derived. Using the definition (3), writing in terms of and |C|'|^, 



using the GOE results \C^\ =1 and \C§\ ln(|C|'| ) = — lnO.48, applying the last equality 
in (1) and replacing all the densities by their corresponding Gaussian forms, converting the 
sum in (3) into an integral and finally carrying out the integration, the expression for Ih in 
wavef unctions is obtained, 



^^^^^ EG OE{l +2) _ J pbiv:giE,Ek) j Pbiv.giE, Ek) 



(6) 



1 — exp — exp — 



. 2 / V 2 ; 

The result in (6) was reported in [|1^] without details. By rewriting the integral in (6) in 
terms of Fk{E) and making small {Sa"^) expansion just as in the case of NPC, the formula 
incorporating corrections due to fluctuations (with respect to k) in the variances of Fk{E) is 
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derived following the arguments that led to (5). Neglecting higher order terms in [(5cr^) / cr|^]^, 



the final result is, 



lHiE) = ^l-Cexpl^ \ exp- ] U 



a 



H 



YiE) \ ; (7) 



Y{E) = {[1- (e' - lY + 4^(1- e) + 2C^ 



IV. NUMERICAL TESTS 



NPC and In are calulated for a EG0E(H-2) with 6 particles in 12 single particle states 
and the results are compared with (4-7) in Fig. 2. In the numerical calculations, the single 
particle energies {i + i = 1,2, . . . , 12 define h{l) (as in 0] and Fig. 1), in the two- 
particle space V{2) is a GOE (calculations use 25 members) with unit matrix elements 
variance and the hamiltonian ensemble is {Ha^x} = ah{l) + A{y(2)}. The value of A = 0.2 
is fixed so that, for a < 1 the level fiuctuations are of GOE; i.e. one is in the chaotic domain 
(see [mplTl and Fig. 1). Results for a = 0, 0.5 and 1 in Figs. 2a,b clearly demonstrate that 
the EG0E(l+2) formulas based on the bivariate Gaussian form for phi^{E, E^) are excellent. 
In these examples C, values are 0.16, 0.59 and 0.82 respectively. The (5cr^) correction is seen 
to be important only for the case with a = 0. In fact the [(5cr^)/o'|/]^ values for the three 
cases considered are 0.121 x 10^^, 0.545 x 10^^ and 0.134 x 10^^. Thus, for realistic fermion 
models that are represented by EG0E(l+2) (with A > A^^.), the correction due to variance 
fiuctuations is expected to be significant only in the situation C, is small. Extension of 
EG0E(2) with explicit inclusion of spin degrees of freedom (each single particle level is 
taken to be doubly degenerate with Sz = ±|; see the third reference in @]) was considered 
and for a system of 6 fermions in 7 levels (i.e. m = 6, = 7 x 2) with total = 0, giving 



d = 1225, NPC was calculated as a function of the excitation energy in [21|; we call this 



model EG0E(2)-S. In this example, as given in [|T|], C = 0.3 and [{Sa^) / crj^] = 0-052. Thus, 
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here the corrections due to variance fluctuations are non-neghgible (the situation in this case 
is similar to the a = case in Fig. 2) and applying (5) gives excellent description, as shown in 
Fig. 3a, of the results for NPC reported in for the EG0E(2)-S model. Returning to Fig. 
2, it should be mentioned that there are differences between the numerical results and the 
predictions based on (4,6) even for the cases with ( = 0.59 and 0.82. These may be due to the 
departures of pbiv{E, Ek) from the bivariate Gaussian form. An important observation from 
(4,6) is, at the spectrum center NPC=((i/3)-\/l — and In = " exp(C^/2). Therefore 
for ('^ close to 0.8 or large, there will be large deviations from GOE even at the spectrum 
center for a system described by EG0E(l+2). This is clearly seen in the a = 1 case in Fig. 
2; here ( = 0.82. Finally, it should be mentioned that the EG0E(l+2) calculations for the 
N = 14, m = 7 system (the case considered in Fig. 1) are also carried out and the results 
are seen to be essentially same as in Figs. 2a,b. 

Let us now turn to the nuclear shell model which is a realistic interacting fermion model. 
There are shell model results for the {2sld) shell (here after called sd shell) nuclei ^^Si 
and ^^Na (see [0] and the second reference in [^) for NPC and Ih in wavef unctions. For 
^^Si the 839 dimensional J = 0,T = space (with six protons and six neutrons in the sd 
shell) and the 3243 dimensional J = 2, T = space are considered. Similarly, for ^^Na the 
307 dimensional J = 2, T = space (with three protons and three neutrons in the sd shell) 
is considered. The results for these nuclei are analyzed using (4,6) as briefly discussed in 
10|,p|. In all the sd shell examples, ( ~ 0.6 — 0.7 and therefore the situation is similar to the 



a = 0.5 case in Fig. 2. Thus, in these examples the departures from GOE at the spectrum 
center are no more than 10% but away from the center, there are large departures. The 
shell model NPC and Ih for sd shell nuclei are seen to be well described by the EGOE forms 
in (4,6). For further confirming this, NPC is evaluated for ^^Mg in the 325 dimensional 
J = 0, T = space (with 4 protons and 4 neutrons in the sd shell) and the results are 
shown in Fig. 3b; here ( = 0.68. It can be concluded that the deviations of the sd shell 
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model results from GOE clearly imply that the shell model hamiltonians are well represented 
by EG0E(l+2) (with A > Xp,.) but not by GOE. It is also seen that the corrections due 
to (5cr^) are small for (sd) shell nuclei (note that here ( is large); in the ^^Mg example, 
[(5cr^)/cr|^]^ = 0.024. In order to further substantiate the EGOE description of the structure 
of nuclear shell model wavef unctions, we have analyzed using (6) the Ih{E) vs E results 
reported recently in [^] for 2plf shell (hereafter called pf shell) nuclei ^°Ca and ^^Sc. In 
the case of ^°Ca the 2051 dimensional J = 6, T = 5 space (with 10 protons in the pf shell) 
and in ^^Sc the 2042 dimensional J = 1,T = 2 space (with one proton and 5 neutrons in 
the pf shell) are considered and a modern large shell model code was used for obtaining the 



Ih values. The shell model results for l}j in Fig. 9 of ||22|, via (2), determine C, to be 0.96 
and 0.92 respectively for the ^°Ca and '^^Sc examples; results for ^^Sc are shown in Fig. 3c. 
From the definition (2) but employing averages over niT spaces (instead of mJT spaces), 
we obtain the C values 0.91 and 0.89 respectively. It should be pointed out that given the 
single particle energies and the two-body matrix elements of the shell model hamiltonians, 
it is easy to calculate C, in fixed niT spaces using trace propagation methods (by extending 
(A. 3) and (A. 4)) ||I],|TB[. The pf shell examples are similar to the a = 1 case in Fig. 2 and 
therefore, as expected, one sees large departures from GOE even at the spectrum center. 
Finally, it is seen from the shell model examples in Fig. 3 and the EGOE examples in 
Fig. 2 that further corrections to the results in (4-7) need to be worked out but this is not 
attempted in this paper. Similarly, study of the nature of fluctuations in NPC and Ih is 
postponed for future. 



V. CONCLUSIONS 

Wavefunction structure given by the EG0E(l+2) random matrix ensemble {H} = 
h{l) + X {V{2)} is studied by deriving compact formulas for NPC and Ih- They are based 
on: (i) the Gaussian form for strength functions Fk{E)^s and the bivariate Gaussian form 
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for pbiv{E,Ek) (with Fk{E) being a conditional density of puv{,E,Ek)) which are vahd in 
the chaotic domain defined by A > Xp^', (ii) there is average- fluctuations separation (with 
little communication between the two) in energy levels and strengths with local strength 
fluctuations following the Porter-Thomas law; (iii) there is a significant unitary group de- 
composition of the hamiltonian. With EG0E(H-2), the NPC and Ih take Gaussian forms 
as a function of the excitation energy and they are defined by the bivariate correlation coef- 
ficient C which measures the variance of the distribution of centroids of Fk{E)'s relative to 
the state density variance. Theory for incorporating corrections due to fluctuations in the 
variances (with k) of Fk{E) is also given. For small (, the present formulation gives back 
the results for pure EG0E(2) (i.e. in the A ^ oo hmit of EG0E(l+2)) as derived in |21[ 
recently. The formulas derived for NPC and Ih are subjected to numerical EG0E(l+2) 
tests with ( changing from 0.1 to 0.8. These and the analysis of the results for a EG0E(2)-S 
example and some nuclear shell model examples, clearly point out that isolated finite real- 
istic interacting particle systems, in the chaotic domain (A > A^?^), will have wavefunction 
structure as given by EG0E(l+2). Finally, in the theory given by (4,6), NPC and Ih depend 
on just one parameter and this appears to be an aspect of 'geometric chaos' (see for a 
recent discussion on 'geometric chaos'). 
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Apendix A 

Let us consider a system of m fermions in N single particle states with a (l+2)-body 
hamiltonian H = h{l) + V{2) where h{l) is specified by the single particle energies ej (with 
i denoting the i^^ single particle state, h{l) = J^i^ini where rij are number operators) and 
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V{2) by the two-body matrix elements Vijki = {kl \ V{2) \ ij). The two-body interaction 
can be seperated into V{2) — -|- V where V^^^ is given by 

V^'^^T^Vijij^nj (Al) 

i<j 

The h{l) + generates the Fk{E) centroids Ek. With N single particle states, there is a 
U{N) group generated by the operators aja^ where aj and Uj are one particle creation 
and destruction operators respectively. With respect to this U (N) group, 1^'°^ decomposes 
into 1/ — 0,1,2 parts and their explicit structure is (for a given m), 

^[01,0 ^ / ^ V ^ r Vo ; 1^0 = 5:%, 

y[0],2 = ^[0] _ y[0],0 _ y[0],l (^_2) 



Similarly the h{l) operator will have u — 0,1 parts; h'^ = me° where e° (A^) ^ and 
^^(1) = X/ ^i^i where e] = ei — e*^. Finally it is to be noted that V behaves as a u = 2 

i 

operator. 

The U{N) norm (in the m-particle spaces) of an operator O is defined by ||C||^ = 



{{O - (C)'")t(C> - (O)™))'". An important theorem is that the i/ = 0, 1, 2 parts of H are 
orthogonal with respect to this U{N) norm. For a, u — 1 operator 0{1) — J^ieiUi, the norm 
square is simply given by 

llfWIll-^^^Eef (A3) 
Similarly ior a u — 2 operator 0{2), 

Using (A. 3) and (A. 4) one can calculate the norms of h} + and and in general 

the later is very small compared to the former. Then h{l) + V^^^ h = J2i Ci^i where 
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= eJ + ^5|Cj (note that at the end we add the spectrum centroid generating part /i^ + y'^l'" 
to h). Thus, neglecting the part, the centroids of Fk{E)^s are generated by h and the 
variances by V. As h and V are orthogonal, afj = + These variances, in m-particle 
spaces, follow easily from (A.3,A.4). See 0,0] for further details. 

Let us consider a EG0E(l+2) hamiltonian H = ah{l) + XV{2) with unit spacing between 
the ej's and the Vijki taken as zero centered Gaussian variables with unit variance. In the 
>> m >> 1 situation one can study the behaviour of ('^ and (5cr^) as follows. The 
correlation coefficient ("^ = and, neglecting the contributions of V{2) to aj^, one 

gets (T^ ~ {mN^/12)a^. Similarly c^y ~ ( ^ ) ( ^ ) A^. Here (A.3) and (A.4) are used. 
Therefore ("^ = [(1 + 3m (A/a)^] ^ and this expression gives 0.51 and 0.76 for the a = 0.5 
and 1 cases in Fig. 2. They compare well with the exact numbers given in Fig. 2. However 
this estimate fails in the situation a — > 0. For a = the h has to be replaced by 

and then the Ek energies are a sum of ^ ™ ^ zero centered Gaussian variables each with 

\ -1 

AT 



variance A . This together with the ay- expression, gives C ~ j j for a ~ as pointed 
out in |21|. The number quoted for the a = case in Fig. 2 is close to this estimate. 
Finally an estimate for is obtained from (A.4) by noting that ay is a sum of 

^ ^ (r)(r) X^"Variables and therefore [(5cr^)/cr|^]^ = 2/K as given first in [|T^. Then, 
a^ = {l- gives the final result [{Sa^)/(rjj]^ ~ 2(1 - 
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FIGURE CAPTIONS 



Fig. 1 Strength functions Fk{E), Dyson-Melita A3 statistic for level fluctuations and oc- 
cupancies {E I rii I E) for EG0E(l+2) for various values of the interaction strength A in 
{H} = h{l) + X{V{2)} for a system of 7 fermions in 14 single particle states (due to compu- 
tational constraints, here only one member is considered just as in |1^); the matrix dimension 
is 3432. The single particle energies used in the calculations are = + = 1, 2, . . . , 14 



just as in [|T^. (a) the histograms are EG0E(l+2) results for the strength functions, con- 
tinuous curves are BW fit and the dotted curves are Gaussian for A < 0.1 and Edgeworth 
corrected Gaussian ||T6| for A > 0.1. In constructing the strength functions, ICf'p are 
summed over the basis states \k) in the energy window E/^ ± A and then the ensemble 
averaged Fji^^{E) vs E is constructed as a histogram; the value of A is chosen to be 0.025 for 
A < 0.1 and beyond this A = 0.1. Here Ek = {E^ — eH)l<^H and in the figure Ey, = 0. Note 
that for ^ 0.2 there is BW to Gaussian transition, (b) The A3(L) statistic for overlap- 
ping intervals of length L < 40 are compared with Poisson and GOE values. For A ~ 0.06, 
there is Poisson to GOE transition in the A3 statistic, (c) The wavy curves are numerical 
EG0E(l-|-2) results for occupancies and the smoothed curves with A > 0.06 correspond to 
the results of EG0E(2) theory (ratio of Gaussians). Note that for A < 0.06 there are wide 
fiuctuations in occupancies and the smoothed forms here are meaningless. All the results 
are shown for the lowest 6 single particle states. Results similar to those in the figure, for 
the A = 12, m = 6 case are reported in [^. 

Fig. 2 (a) Number of principal components NPC and (b) the localization length In in 
wavefunctions for a system of 6 interacting particles in 12 single particle states (matrix 
dimension is 924). Here, for conveniance, the EG0E(l+2) hamiltonian is changed to 
{H(^a,x)} = ah{l) + A{y(2)}. Numerical EG0E(l+2) results correspond to filled circles. 
The continuous curves correspond to the theory (4) for NPC and (6) for Ih- For the case 
with a = 0, the dashed curves correspond to the theory (5) for NPC and (7) for Ih- For 
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other cases the correction due to variance fluctuations is neghgible and hence only the results 
of (4,6) are shown in the figure. Note that NPC=(i/3 and = 1 for GOE. See text for 
further details. 

Fig. 3 (a) NPC for the EG0E(2)-S model described in the text compared with the results 
given by (4,5). The filled circles are for the numerical EG0E(2)-S calculations reported in 
The continuous and dashed curves represent (4) and (5) respectively, (b) NPC for the 
{sd) shell nucleus ^^Mg compared with (4). The shell model calculations are same as in . 
In this example (4) and (5) give almost identical results and hence the curve corresponding 
to (5) is not shown in the figure, (c) Shell model results for ^^Sc for Ih reported in ||2^ 
compared with the theoretical curve given by (6) with = 0.92. 
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